Two-dimensional Bose gas under extreme rotation 
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We show that a Bose-condensed gas under extreme rotation in a 2D anisotropic trap, forms a 
novel elongated quantum fluid which has a roton-maxon excitation spectrum. For a sufficiently 
large interaction strength, the roton energy reaches zero and the system undergoes a second order 
quantum transition to the state with a periodic structure - rows of vortices. The number of rows 
increases with the interaction, and the vortices eventually form a triangular Abrikosov lattice. 
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Rotating Bose Einstein condensates (BEC) of trapped 
atoms constitute a novel many-body system where nucle- 
ated quantized vortices form a triangular lattice Q , and 
a fast rotation is expected to change dramatically the 
properties of the gas. In a harmonically trapped BEC 
rotating at a frequency £1 close to the trap frequency in 
the rotation plane u±, the vortex lattice can melt when 
the number of vortices approaches the number of parti- 
cles 0,13- In this respect, there is an analogy with type-II 
superconductors, which undergo a transition to the nor- 
mal phase above a critical value of the magnetic field 
4] . Possible signatures of melting of the vortex lattice in 
rapidly rotating BEC's were observed in experiments at 
JILA and at ENS [3 . The proposals that are put forward 
to describe the state of the rapidly rotating Bose gas in- 
clude yrast states 0, correlated quantum Hall states of 
bosons 0, , and a giant vortex state [j| . 

In this letter we show that an anisotropy in the trap- 
ping potential can drastically change the picture. As 
found at ENS Q, a Bose gas under critical rotation 
(Q ~ lui_) can become very elongated in one direction 
in the rotating frame. We focus on this case and assume 
that f£ is close to the smallest of the trap frequencies 
in the rotation plane. Then the condensate becomes free 
along the direction of the weaker confinement and forms a 
novel quantum fluid in a narrow channel. The excitation 
spectrum of this fluid has a " roton-maxon" character and 
becomes unstable at a critical interaction strength. This 
instability leads to the formation of a periodic structure 
which represents vortex rows. An increase in the inter- 
action (decrease in the anisotropy) increases the number 
of rows and reduces the correlation length, and the gas 
ultimately enters the strongly correlated regime. 

We consider a two-dimensional (2D) Bose-condensed 
gas at zero temperature, rotating with frequency £1 and 
harmonically trapped with frequencies uj XtV — oj± \f \ =F e 
along the x,y axes in the rotating frame. The Hamilto- 
nian of the system in this frame reads (see [l(ij]): 



H= d 2 r¥ 



(-ihV+mrxQ) 2 
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where W ,Si! are bosonic field operators, m is the parti- 
cle mass, g is the coupling constant for the mean-field 
interaction [llj, and the effective trapping potential is: 



V ef f(r) = (m/2)[(f4 - Vl 2 )x 2 + (cj 2 y - tf)y 2 



(2) 



The Hamiltonian JJ) is analogous to that of charged par- 
ticles in the magnetic field, and in this respect the quan- 
tity mf x Cl is the gauge field. We consider the limit of 
extreme rotation, where f2 = oj x = wxvl — e and the gas 
becomes free along the x direction. Then, assuming that 
in this direction the atoms are confined in a large rectan- 
gular box of size L, the gas becomes a long cigar. After 
the gauge transformation ^ = •q/ e lmnx v/ h ^ the Hamilto- 
nian in the Landau gauge can be written as: 



H= d 2 r& 



( Px + 2mny) 2 +p 2 
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where cj_ = cjj_v2e -C fl, assuming a small ellipticity e. 

Omitting the interaction term, the single particle 
eigenstates are the Landau levels [l2J separated from each 
other by an energy gap ~ 2MI. In the dilute limit, where 
the mean field interaction gri2D Ml, {n^D is the 2D 
density), we may restrict our discussion within the lowest 
Landau level. Then the field operator can be written in 
the form ^ = 4>k^ki where ctfc is the creation operator 
of a particle with momentum k along the x direction and 
4>k is the corresponding eigenfunction: 



<t>k{x,y) 



exp(zfcx) 



exp 



(4) 



where fi = ^Jtt 2 +u 2 _/A, and l = (h/2mQ) 1 / 2 . Then, 

after the spatial integration of Eq. we obtain an ef- 
fective one-dimensional (ID) Hamiltonian: 



H 



{h 2 k 2 /2m*)4a fc + {g*/2L) J2 4 

k,k' ,q 

exp{-l 2 Q [(k-k' + q) 2 +q 2 }/2}, 
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where g* = g/^/2~rrlo is an effective ID coupling constant. 
The Hamiltonian © describes particles with a large ef- 
fective mass m* = m(2il/w_) 2 ^> m. The fact that m* 
is not infinite and the kinetic energy term is still present 
originates from the asymmetry of the trapping potential. 
This asymmetry leads to a small difference between the 
frequencies Q and f2. However, once the finite kinetic en- 
ergy term is extracted one may put Q = Q, and we have 
done this in the second term on the rhs of Eq. © . 

The momentum dependence of the interaction term 
in the Hamiltonian © originates from the presence of 
the gauge field. The wavefunctions of particles which 
have opposite momenta in the x direction are shifted in 
opposite directions along the y axis, which decreases their 
overlap and reduces the interaction amplitude. 

The behavior of the system is governed by the particles 
with momenta k < Iq 1 , for which the extension of the 
wave function in the y direction is Iq. If the ID density 
n = N/L (N is the total number of particles) satisfies 
the condition uIq -C 1, then we are dealing with a ID 
Bose gas. In this case, characteristic momenta that are 
important are at least of the order of n^ 1 or smaller. 
They satisfy the condition Mq *C 1, and the exponential 
term in Eq.© is equal to unity. Then the Hamiltonian © 
corresponds to the Lieb-Liniger model for the ID Bose 
gas, well described in literature [ljj. 

We, therefore, focus on the other extreme, where 

nl > 1. (6) 

Then the system can be viewed as a 2D gas in a narrow 
channel. The 2D density is ~ n/lo and the interaction 
energy per particle is I ~ ng/lo. The characteristic ki- 
netic energy of a particle at the mean distance from other 
particles is K ~ h 2 n/m*lo, and for K 3> / the wave func- 
tion of the particle at such interparticle distances is not 
influenced by the interactions. The gas is then in the 
weakly interacting mean-field regime. The criterion of 
weak interactions takes the form (see, for example [lij): 

m*g/h 2 < 1. (7) 

The correlation length is l c ~ h/^/m*g*n, and for l c 3> Iq 
the gas enters the ID regime. In this case the ID criterion 
of weak interactions, m*g* /h 2 n <C 1, is automatically 
satisfied under the conditions © an d ©• 

For the weakly interacting 2D Bose gas in a narrow 
channel, as well as in the ID Bose gas, for a sufficiently 
large size L the ground state can be a quasi-condensate. 
In this state the density fluctuations are suppressed, but 
the phase fluctuates in the x direction on a distance 
scale greatly exceeding l c Q. However, locally the 
quasi-condensate is indistinguishable from a true BEC. 
As follows from the analysis of hydrodynamic equations 
in the density-phase representation [14j . the excitation 
spectrum is the same as the one obtained in the Bogoli- 
ubov approach assuming that most particles are in the 



condensate. Employing this approach we first reduce the 
Hamiltonian © to a bilinear form: 

H b = [h 2 k 2 /2m* + 2n 9 *exp{-fc 2 Z 2 ,/2}] a{a k 

k 

+ (ng*/2)Y,exp{-k 2 l 2 }(ala f _ k +a k a^ k ) (8) 

k 

Then, diagonalizing the Hamiltonian © we obtain the 
excitation spectrum: 

e 2 (k) = [h 2 k 2 /2m* +g*n{2e^p(~k 2 l 2 /2) - l}f 

-.g* 2 n 2 exp{-2fc 2 /^} . (9) 

The key feature of the excitation energy 10 is the 
momentum dependence of the interaction terms propor- 
tional to g* . The structure of the spectrum depends on 
the ratio of the mean-field interaction to the kinetic en- 
ergy at momentum k w 1/Iq. This ratio can take both 
small and large values and is given by: 

In units of mMl/m*, the excitation energy is a univer- 
sal function of (3 and kl$. For small (3, the interaction 
terms proportional to g* in Eq.© are important only 
at k <C 1/io; where they become momentum indepen- 
dent. Then Eq.© gives the ordinary Bogoliubov spec- 
trum, with a small sound velocity c s = y/ng* /m* . 
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FIG. 1: Excitation energy (in units of mMl/m*) versus Mo. 

For P > 1 the situation drastically changes. The in- 
teraction is already important for momenta k > 1/Iq, 
where the interaction dependent terms in Eq. © decrease 
with increasing k. For this reason the spectrum develops 
the roton-maxon structure for (3 > 2.6 (see Fig.l). This 
structure is well-known in the physics of superfluid 4 He 
and is predicted for trapped dipolar condensates |l5l | and 
studied for lattice bosons with long-range interactions 
|l6| . The roton minimum is located at k ~ 1/Iq and the 
corresponding excitation energy decreases with increas- 
ing (3. For a critical value (3 = 4.9, the roton minimum 
reaches zero at k = k c = 1.6/Zo aud a further increase in 
(3 makes the system unstable. We thus conclude that our 
weakly interacting 2D Bose-condensed gas in a narrow 
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channel is stable for (3 < 4.9, exhibiting a roton-maxon 
spectrum in the range 2.6 < (3 < 4.9. For (3 > 4.9, where 
the Bose-condensed state is unstable, one has to find a 
new ground state. 

At the instability point the excitations with momenta 
±fc c can be excited without any cost of energy. This 
indicates that the ground state macroscopic wavefunction 
can contain several momentum components. A general 
form of this type of wavefunction reads: 



FIG. 2: The number of components (N r ) in the ground state 
macroscopic wavefunction for various values of (3. 



ip = VN 



C <Po + X](CM fci + C-kA- 



(11) 



The absence of current in the x-direction requires | C^ | = 
|C_fcJ, and the normalization condition reads |Co| 2 + 

2 Ei ICfeJ 2 = !• Notc tnat Eq.(HIl) gives two possibili- 
ties: a (2j + 1) component wavefunction with Co =/= 0, 
and a 2j component wavefunction for which Cq — 0. 

To understand the instability of the single component 
state, we consider a wavefunction with three components: 



(12) 



$ = VN[C k <p k + C o e l& o + C_ fc 0_ fe ] 



Here all the C-coefficients are real, and we may put 
Cfe = C-k- The quantity \Ck\ 2 represents the occupation 
number of the mode with momentum k. Normalization 
of the wavefunction requires 2|Cfc| 2 + C 2 — 1. We then 
find the critical value (3 C of the parameter (3, above which 
this 3-component state has lower energy than the single 
component one. The wavefunction (|12fl leads to eto = Co, 
a±k = Ck in the HamiltonianJSJ), and we obtain the en- 
ergy per particle: 

E[C k , k]/N = A(k)\C k \ 4 + B(k)\C k \ 2 + g*n/2, (13) 

where the coefficients A(k) and B(k) are given by: 

A(k) = g*n[3-m + 2t 4 -4t 2 cos(29)} (14) 
B(k) = h 2 k 2 /m* -g*n[2-4t- 2£ 2 cos(26% (15) 

with t = exp{-fc 2 Z 2 /2}. For |C fc | 2 < 1/2, the energy 
is minimized for cos(2#) = —1, and we have A(k) > 
0. Minimizing E with respect to | C fe | 2 , yields | C fc | 2 = 
— B(k)/2A(k). The energy of the 3-component state is 

E/N = g*n/2- B 2 /4A, 

and it is lower than the energy of the single compo- 
nent state. The physically acceptable solution requires 
B(k) < 0. Therefore, the transition from the single to 3- 
component state occurs when the minimum value of B(k) 
reaches zero. Using Ea. (ll5|) we find that this happens at 
k = k c and (3 = (3 C = 4.9. So, the 3-component solution 
becomes the ground state for f3 > 4.9. This breaks the 
translational symmetry and gives rise to a modulation of 
the density along the cc-axis with a period of 27r/fc c , which 



is similar to a scenario proposed for superfluid 4 He flow- 
ing along a capillary with a velocity exceeding the critial 
Landau velocity |l7l |. 

A 3-component wavefunction can be viewed as two vor- 
tex rows along the x-axis. The nodes in the x, y plane are 
obtained straightforwardly, and near the transition point 
they are very far from the line y = 0. The transition from 
the single to 3-component wavefunction can be treated 
as a second order quantum transition. The energy and 
chemical potential [i change continuously, whereas the 
compressibility undergoes a jump at the transition point. 
For the 3-component state it is smaller by an amount 



A['^| = 

on 



2.9* [1 - t c 



[3 - 8t c + Ml + 2ti) 



(16) 



where t c = exp{-fc 2 Z 2 /2} w 0.28. 

Near the transition point quantum fluctuations in- 
crease due to the vanishing excitation energy at the roton 
minimum. This energy can be expressed as 



e{k) 



2m*ll 



(f3 c -f3) +1 lt{k 2 -k 2 c ) 2 

Pc 



n 1/2 



where 7 w 0.12. For the mean square fluctuations of the 
density we then obtain: 



dk 
2mr 



h 2 k 2 



2m*e(k) 



1 



1 



nlo \Pc—/3 



These fluctuations become large for 8(3 — ((3 C — f3) < 
exp(— nlo). Thus, the transition from a single to 3- 
component wavefunction occurs in the interval 8(3, which 
is exponentially narrow due to the inequality ©. The 
behavior of the system in this region requires a special 
investigation and is beyond the scope of this paper. 

It is important that already for (3 — 5.4 the ground 
state wavefunction changes from 3 to 2-component, and it 
again becomes 3-component at (3 = 20. With increasing 
(3 (either increasing g*n or decreasing more momen- 
tum states are macroscopically occupied. This is because 
an increase of (3 is equivalent to increasing the effective 
mass to*, which makes momentum states in the lowest 
Landau level (LLL) more degenerate. 

Our results for the number of momentum components 
in the ground state wavefunction are displayed in Fig. 2. 
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FIG. 3: Density profile of 2-component (I), 3-component (II), 
and 5-component (III) states. (Lengths are in units of Iq). 

These states describe one or several vortex rows along 
the x-axis. For example, a two component state repre- 
sents one vortex row (see Fi g. 3 I). Eventually, a triangu- 
lar Abrikosov vortex lattice |l8j is formed when increas- 
ing (3 to very large values. It should be mentioned here 
that the structure of vortex rows has been obtained in the 
studies of type-II superconductors an d m the studies 
of condensates in rotating anisotropic traps |2(i| . 

In the presence of many rows of vortices, we may 
rely on an average vortex description in the LLL. Us- 
ing the relation < 4>k\y 2 / lo\<l>k >= k 2 l$ + 1/2, we write 
the energy functional of the system in the form: E[p] — 
J d 2 r [muj 2 _y 2 p + gp 21 \ /2. Minimizing the energy we ob- 
tain the coarse grained density p(y) — ^ [/i — ^mu} 2 _y 2 ^ , 
which smooths out density modulations introduced by 
vortex rows. The size of the vortex core is always ~ lo 
and the number of the vortex rows is ~ /3 1 / 3 . The energy 
per particle is given by: 

f^(s)l< 3 ^/") 2/3 ' <"> 

The result of Eq. (|17fl deviates by less than 10% from the 
energy obtained by the full numerical minimization. 

The vortex lattice is expected to melt when the num- 
ber of vortices N v approaches the number of particles N 
In our case this condition leads to m*g/h 2 ~ n 2 l 2 . 
However, our approach of the weakly interacting gas re- 
quires m*g/H 2 -C 1 and, as we are considering large val- 
ues of n/o, it breaks down before the melting transition. 
Therefore, the description of this transition requires a 
more elaborate treatment of quantum fluctuations. 

In conclusion, we have shown that a 2D BEC at ex- 
treme rotation frequency in an elliptic trap forms a novel 
quantum fluid in a narrow channel. The behavior of this 
fluid is determined by the parameter Eq. (|10fl which in- 
creases with the interaction strength and effective mass 
m* . For (3 < f3 c — 4.9, the excitation spectrum of the 
fluid has a small sound velocity and exhibits a roton- 
maxon character. At a critical interaction strength /3 C , 
the roton energy reaches zero and the uniform (in the 
long direction) ground state becomes unstable. The sys- 
tem then undergoes a second order quantum transition to 
the state with a periodic structure which can be viewed 



as rows of vortices. For f3 > j3 c , with increasing the in- 
teraction parameter /3, more vortex rows are nucleated. 
Finally the vortices form the Abrikosov lattice which can 
melt due to quantum fluctuations. 
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